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Microscopic theory of Gilbert damping in metallic ferromagnets 
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We present a microscopic theory for magnetization relaxation in metallic ferromagnets of 
nanoscopic dimensions that is based on the dynamic spin response matrix in the presence of spin- 
orbit coupling. Our approach allows the calculation of the spin excitation damping rate even for 
perfectly crystalline systems, where existing microscopic approaches fail. We demonstrate that the 
relaxation properties are not completely determined by the transverse susceptibility alone, and that 
the damping rate has a non-negligible frequency dependence in experimentally relevant situations. 
Our results indicate that the standard Landau-Lifshitz-Gilbert phenomenology is not always ap¬ 
propriate to describe spin dynamics of metallic nanostructure in the presence of strong spin-orbit 
coupling. 


Magnetization relaxation in metals is at the heart of 
spin current generation and detection processes currently 
under investigation, many of them candidates to play 
protagonist roles in innovative spintronic devices. The 
Landau-Lifshitz-Gilbert (LLG) equation is widely used 
to describe the spin dynamic properties of magnetic ma¬ 
terials mm- It includes an important system-dependent 
parameter, called the Gilbert damping constant, usually 
denoted by ac, that regulates the relaxation of the mag¬ 
netization towards stability, after it is driven out of equi¬ 
librium. Recently, a lot of effort has been put into the 
determination of this damping rate Hi], which charac¬ 
terizes the pumping and absorption of pure spin currents 
in nanostructures that are of great interest in the field of 
spintronic. In most of them spin-orbit interaction is sig- 
nihcant, and responsible for a desirable interplay between 
charge spin and angular momentum excitations. 

There is a general agreement between practitioners in 
the field that a proper microscopic theory of magnetiza¬ 
tion relaxation in metals requires a good description of 
the electronic structure of the system including spin- 
orbit coupling [3H8]. The conventional approach is to 
combine a realistic electronic structure with some kind of 
adiabatic approximation to derive expressions that can 
be directly related to the Landau-Lifshitz-Gilbert phe¬ 
nomenology. This strategy has been employed by Kam- 
bersky [3] and many others since [MS]. This conven¬ 
tional approach has important limitations. It neglects the 
coupling between transverse spin, longitudinal spin and 
charge excitations (which is an important consequence of 
the spin-orbit coupling), and incorrectly pedicts the di¬ 
vergence of the damping parameter for a perfectly crys¬ 
talline system. Actually, for ferromagnets that display 
rotation symmetry in spin space, the Goldstone theorem 
ensures that any experiment which measures the total 
transverse magnetic moment of the sample will produce 
a resonant response with zero linewidth |9] . In the pres¬ 
ence of spin-orbit interaction, however, this symmetry is 
explicitly broken, and the resonant spectrum acquires a 
finite linewidth m- 

We put forward a more fundamental microscopic ap¬ 
proach to the calculation of the spin dynamics damp¬ 


ing rate that takes fully into account the effects of SOG 
on the spectrum of spin excitations of itinerant sys¬ 
tems. Namely, we consider the coupling of transverse 
spin excitations to longitudinal spin and charge excita¬ 
tions, induced by the spin-orbit interaction. We calculate 
the FMR spectrum at finite frequencies and arbitrary 
anisotropy values, without employing any adiabatic ap¬ 
proximation. We will show that those ingredients are es¬ 
sential to correctly describe the magnetization relaxation 
in very clean metallic ferromagnets of nanoscopic dimen¬ 
sions, and that the Landau-Lifshitz-Gilbert phenomenol¬ 
ogy fails to capture essential features of the magnetiza¬ 
tion dynamics in those systems. 

This letter is organized as follows: we will present 
briefly our formalism, discuss its main features and 
present numerical results for two model systems that il¬ 
lustrate common but qualitatively different situations. 

General Formalism - The spectrum of spin excitations 
of a ferromagnet can be obtained from the spectral den¬ 
sity associated with the transverse spin susceptibility 

X+-(Z,r;II)= |dte“((5+(t),5,7(0))), (I) 

where 

{{s+{t),s^m) ^ -z9ms+{t),si;m), ( 2 ) 

and 

~ X (^) 

n 

The operator creates one electron in the atomic ba¬ 
sis state ^ localized at lattice site I with spin a. Although 
we are usually interested in x~^ ihl': as defined above, 
its equation of motion involves the orbital-resolved sus¬ 
ceptibility. 

In the absence of spin-orbit coupling (SOG) and within 
the random phase approximation (RPA), the equation of 
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motion for closed and x"*" {l,l'-,V,) can 

be expressed in the well-known RPA form, 

x+-(r!) = [1 + ( 5 ) 

where Xo (^) is the mean-field (sometimes called non¬ 
interacting, or Hartree-Fock) susceptibility. This expres¬ 
sion is schematic and must be understood as a matrix in 
orbital and site indices, in real space, or a wave-vector 
dependent matrix in reciprocal space. The crucial point, 
however, is that, in the absence of spin-orbit coupling, 
within the RPA, the transverse spin susceptibility is un¬ 
coupled from any other susceptibility. This ceases to 
be true when SOC is included, as we demonstrated in 
ref. [ini x^ becomes coupled to three other susceptibil¬ 
ities, namely 



(7) 
■ ( 8 ) 

The system of equations of motion obeyed by these 
four susceptibilities can be cast into a form strongly re¬ 
sembling the RPA result by introducing a block-vector 
X = with x^^^ = With 

an equivalent definition for the mean-field susceptibili¬ 
ties Xo™^ we write 

x(f^)=Xo(^^)-Ax(F!), (9) 

where the “super-matrix” A is proportional to the effec¬ 
tive Coulomb interaction strength and involves convolu¬ 
tions of single particle Green functions. Explicit forms 
for its matrix elements are found in Ref. nUl The nu¬ 
merical analysis of the susceptibilities x^^K X^^^ ^ind x^^^ 
show that their absolute values are many orders of mag¬ 
nitude smaller than those of x^^^ = X"' ■ It is, thus, 

tempting to argue that the transverse susceptibility is ap¬ 
proximately decoupled from x^‘^\ X^^^ x^^^ and that 

it can be calculated via the usual RPA expression with 
the single particle Green functions obtained with spin- 
orbit coupling taken into account. This is not a good 
approximation in general, since the matrix elements of A 
that couple x^^^ to the other susceptibilities are far from 
negligible. Our numerical calculations indicate that they 
are essential to determine correctly the features of the 
FMR mode around the resonance frequency. Thus, the 
behaviour of x^^^ in the presence of spin-orbit coupling 
cannot be inferred from Xo^^ in the zero-frequency limit 
alone, as it is usually assumed in the literature on the cal¬ 
culation of the Gilbert damping parameter msiaiii]. 

Numerical Results - We start the discussion by pre¬ 
senting results for the Gilbert constant ac for unsup¬ 
ported ultrathin Co films. Here we determine ao from 
the ratio between the FMR linewidth Ail and the reso¬ 
nance frequency Hq. First we turn off spin-orbit coupling 


FIG. 1: Gilbert damping constant ao as a function of the 
imaginary part rj added to the real energy, for an ultrathin 
film of two atomic layers of Co where SOC has been turned 
off. It is clear that ao vanishes as 77 —>■ 0. 

to check the consistency of our approach. Even with 
SOC turned off we still find a finite linewidth for the 
EMR mode. It comes, as we will shortly demonstrate, 
from the small imaginary part rj that is usually added 
to the energy in the numerical calculations of the sin¬ 
gle particle Green functions, in order to move their poles 
from the real axis. We calculate ao for various values of 
r] and extrapolate to 77 —>■ 0+, as shown in Fig. It is 
clear that lim^_,,Q+ ac = 0. Thus, our approach correctly 
predicts that the Gilbert damping constant vanishes in 
the absence of SOC, as it should. Indeed, it is easy to 
show [3] that the FMR mode is a stationary state of the 
mean-held hamiltonian and, as such, has inhnite lifetime 
in the limit 77 —> 0^. Now we discuss the dependence 
of ac on rj for a hxed, non-zero value of the spin-orbit 
coupling strength We used LG AO parameters appro¬ 
priate for bulk Co to describe the electronic structure 
of all Co hlms we investigated. The quantitative details 
of the ferromagnetic ground state and excitation spec¬ 
tra are sensitive to the LCAO parameters used, but their 
qualitative behaviour is very robust to small changes in 
the electronic structure. Our strategy is to use the same 
set of LCAO parameters for all him thicknesses to avoid 
modihcations in ao coming directly from changes in the 
LCAO parameters. This allows us to focus on geometric 
effects and on the 77 -dependence. 

Figure [^shows the dependence of the Gilbert damping 
constant ac on the imaginary part 77 for Co hlms of var¬ 
ious thicknesses. Clearly ac approaches hnite values as 
77 —0. Cobalt has a small spin-orbit coupling constant. 
We would like to investigate the effect of increasing the 
strength of the SOC on the damping rate. Instead of 
artihcially increasing ^ in Co we consider a more realis¬ 
tic setting where a double layer of Co is attached to a 
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FIG. 2: Gilbert damping constant ac as a function of the 
imaginary part 77 added to the energy, for Go ultra thin films 
of various thicknesses: 1 (circles), 2 (squares), 4 (diamonds) 
and 6 (triangles) atomic layers. The strength of the SOG is 
^ = 85 meV. The solid lines are guides to the eye. 


non-magnetic substrate with high SOC parameter, such 
as Pt. This system has a particularly interesting fea¬ 
ture: the magnetization easy axis is perpendicular to the 
plane. However, we found that, for the LCAO param¬ 
eters we employed, the magnetization in-plane is also 
a stable configuration, with a small magnetocrystalline 
anisotropy. The damping rate, however, is much larger 
in the 2Co/2Pt system than in the unsupported Co films. 
This is a nice example of how the anisotropy energy is 
strongly influenced by the system’s symmetry, but the 
damping rate is relatively insensitive to it, depending 
strongly on the intensity of the spin-orbit coupling. It is 
also an extremely convenient situation to test an assump¬ 
tion very frequently found in the literature on Gilbert 
damping, although sometimes not explicitly stated: that 
the FMR linewidth Ail is linearly dependent on the res¬ 
onance frequency Hq and that AH —0 as Oq —>■ 0. This 
is not an unreasonable hypothesis, considering the weak 
static helds commonly used in FMR experiments and the 
smallness of the spin-orbit coupling constant, compared 
to other energy scales of a ferromagnet. Our calculations 
for the Co films confirm that this relationship is approx¬ 
imately held. In this case, the Gilbert constant aa may 
be extracted from the FMR spectrum by simply fitting 
it to a Lorentzian and is practically field-independent. 
However, our results for 2Co/2Pt indicate that the FMR 
linewidth is finite as Hq —> 0 , leading to a significantly 
frequency-dependent ao, as shown in Fig. In order to 
illustrate how the determination of a damping parame¬ 
ter is affected by the finite value of AH as Hq —>■ 0 we 
extracted the linewidths from the calculated spectra for 
the 2Co/2Pt system by fitting Lorentzians to our calcu¬ 
lated spectral densities. The results are shown in Fig. 
One of its most important consequences is that, if one 
wishes to define a value of ao for the system above, it 
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FIG. 3: a) Spectral densities of the FMR mode for the 
2Go/2Pt system subjected to various static magnetic fields 
(from -0.3 T to 4 T). The inset shows the resonance frequency 
as a function of the Zeeman field B. b) The Gilbert damping 
parameter aa as a function of applied Zeeman field B. The 
inset shows the FMR line width as a function of resonance 
frequency Hq. The strengths of the SOC are ^co = 85 meV 
and ^pt = 600 meV. 


must be defined as a function of the Zeeman field, as is 
illustrated in Fig. In principle this poses a problem 
for the procedure usually employed to determine FMR 
spectra experimentally, since there the free variable is 
the Zeeman field, not the frequency of the exciting field. 
In Fig. 1^ we illustrate this issue by plotting the FMR 
spectral density as a function of the Zeeman held for two 
hxed pumping frequencies, 24 GHz and 54 GHz. The 
curves have nice Lorentzian shapes, but the values for 
the Gilbert damping parameter ao extracted from these 
curves depend on the pumping frequency {ao = 0.034 
for Hq = 0.10 meV and ao = 0.042 for Hg = 0.22 meV). 
Also, they do not correspond to any of the values shown 
in Fig. io, although the Zeeman held values that de¬ 
termine the linewidth in Fig. lie within the range of 
Zeeman held values showed in Fig. HJd. Thus, if ao is 
dehned as AH/Hq, its value for a given sample depends 
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FIG. 4: Spectral densitty of the FMR mode for the 2Co/2Pt 
system plotted as a function of the Zeeman field B at a 
fixed pumping frequencies: Vp = 24 GHz (squares) and 
Vp = 54 GHz (circles). The solid curves are Lorentzian fits to 
the calculated points. 
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FIG. 5: a) Spectral densities of the FMR mode for the 
2Go/2Pt system with Stoner enhancement in Pt turned on 
(black line) and off (red line). 


on wether the FMR spectrum is obtained in a fixed fre¬ 
quency or fixed Zeeman field set ups. Our results also 
imply that the existing expressions for the damping con¬ 
stant ac are not valid in general, specially for very clean 
systems with large spin-orbit coupling materials. The 
conventional approaches express ac as the ratio AO/Oq 
in the Oq —^ 0 limit. As we have just shown, this limit 
does not exist in general, since AO approaches a finite 
value as Oq —0 . 

In experimental papers [12l HI] the FMR linewidth is 
assumed to have a zero-frequency offset, just as we de¬ 
scribed. This is usually attributed to extrinsic broad¬ 
ening mechanisms, such as two-magnon scattering El: 
due to the combination between inhomogeneities in the 
magnetic films and dipolar interactions. This is certainly 
the case in systems with small SOC, such as Fe films de¬ 
posited on GaAs or Au m- However, we have shown 
that there can be zero-frequency offset of intrinsic ori¬ 
gin if the SOC is large. The effect of this intrinsic offset 
should be easily separated from that of the two-magnon 
scattering mechanism, since the latter is not active when 
the magnetization is perpendicular to the plane of the 
film [T^ . 

We would like to remark that Stoner enhancement in 
Pt plays a very important role in the determination of the 
damping rate. We had shown previously m that, in the 
absence of spin-orbit coupling, Stoner enhancement had 
a very mild effect on the damping rate in the Co/Pd(00I) 
system. In the presence of SOC, however, the effect can 
be very large indeed. Both magnetocrystalline anisotropy 
and damping rate are significantly different in the en¬ 
hanced and non-enhanced cases, as shown in Fig. The 
Gilbert parameter is also very different in the two cases: 
a™*' = 0.11, whereas = 0.33. Thus, proper 


treatment of Stoner enhancement in substrates like Pd 
an Pt is essential for the correct determination of spin 
relaxation features. 

We presented a microscopic approach to the calcu¬ 
lation of the Gilbert damping parameter ac for ultra- 
thin metallic magnetic films, illustrated by results for 
Co films and Co/Pt bilayers. Our approach is based on 
the evaluation of the dynamic transverse susceptibility in 
the presence of spin-orbit coupling, taking into account 
realistic electronic structures and the coupling between 
transverse spin, longitudinal spin and charge excitations. 
It predicts finite values of ac in the limit of perfectly 
crystalline films, a regime where methods based on the 
torque correlation formula find a diverging Gilbert damp¬ 
ing parameter. We showed that the coupling between 
transverse, longitudinal and charge excitations, due to 
spin-orbit coupling, is of fundamental importance for the 
correct determination of FMR spectra in metallic sys¬ 
tems. We have also shown that the damping rate ex¬ 
tracted from the FMR spectrum for fixed pumping fre¬ 
quency differs considerably from that extracted from the 
FMR spectrum for fixed Zeeman field. In this case the 
Gilbert damping parameter ac becomes frequency de¬ 
pendent, in contrast to what is assumed in the standard 
Landau-Lifshitz-Gilbert phenomenology. Moreover, we 
have numerical indications that the Gilbert parameter 
is not well defined in the limit of vanishing resonance 
frequency, a fact that is very relevant to calculational 
schemes based on the adiabatic approximation. Inci¬ 
dentally, Stoner enhancement in materials like Pt and 
Pd also plays an important role in the determination of 
FMR frequencies and damping rates. These results may 
lead to important modifications of the interpretation of 
damping “constants”, either calculated or inferred from 
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experimental results, for systems where spin-orbit cou¬ 
pling is strong. We believe these issues may be crucial 
for the correct description of relaxation in very clean sys¬ 
tems of nanoscopic dimensions, specially in the presence 
of relatively weak magnetocrystalline anisotropy. 
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